Ground state of finite nuclei evaluated from realistic interactions 
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Ground state properties of finite nuclei ( 16 and 40 Ca) are evaluated from 
realistic nucleon-nucleon interactions. The calculations are based on the 
Brueckner-Hartree-Fock approximation. Special attention is paid to the role 
of the energy spectrum for the particle states, in particular for those close to 
the Fermi energy. Additional binding energy is obtained from the inclusion 



■ of the hole-hole scattering term within the framework of the Green function 



approach. Results for the energy distribution of the single-particle strength 
and the sensitivity to the nucleon-nucleon interaction are investigated. For 



C$ [ that purpose three modern nucleon-nucleon interactions are employed: the 

Argonne V18, the charge dependent Bonn potential and a realistic nucleon- 
qq ■ nucleon interaction which is based on chiral perturbation theory and which 

has recently been fitted by the Idaho group. 
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I. INTRODUCTION 



One of the central aims of theoretical nuclear physics is the attempt to determine the bulk 
properties of nuclear systems like its binding energy and its density or radius from a realistic 
model of the nucleon-nucleon (NN) interaction, i.e. a NN interaction which yields an accurate 
fit of the NN scattering data below the threshold of 7r-production. Quite some progress has 
been made during the last ten years in the determination of such realistic NN interaction 
models. A family of charge- dependent NN potentials has been generated, which all fit the 
empirical scattering data with high precision 0-0] • Nevertheless, using such NN interactions 
in nuclear structure calculations, one obtains results, which exhibit significant differences 



The NN interactions containing non-local terms like the charge-dependent Bonn potential 
0, CDBonn, or one of the Nijmegen interaction model |1[ tend to be "softer" than the purely 
local interactions like the Argonne V18 potential @. Here the expression softer interaction 
is used to identify those interactions which induce weaker two-nucleon correlations in the 
nuclear many-body wave function. This implies that the total energy of nuclear matter 
calculated at a given density without correlations, i.e. using the Hartree-Fock approxima- 
tion, yields less repulsive results for a soft, non-local potential as compared to the energy 
calculated for a stiffer, local NN interaction. As an example we mention that Hartree-Fock 
calculations for nuclear matter at the empirical saturation density yield 4.6 MeV per nucleon 
using the CDBonn interaction while 30.3 MeV are obtained if the Argonne V18 potential is 
used ||. Including effects of NN correlation as it is done e.g. in a Brueckner-Hartree-Fock 
(BHF) calculation leads to results which are attractive and much closer to each other. 
All these modern interactions, however, tend to predict too much binding energy for nuclear 
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matter and saturation densities which are too large. This is the case for the BHF approach 
but is also observed if variational calculations are performed |],[7j Softer non-local interaction 
yield larger binding energies and saturation densities than local interactions ||. This over- 
binding is often compensated by a three-nucleon force, which can be adjusted in such a 
way that the empirical saturation point of nuclear matter is reproduced []8| p~2[] . These three- 
nucleon forces can be understood to simulate the relativistic effects as incorporated e.g. in the 
Dirac-BHF approach [|X5| P~5|] . These three-nucleon forces, however, can also be interpreted 
to describe the effects of virtual excitations of nucleons to the A(3,3) [fL7| , |l8[| or TV* Roper 
resonance P,|H?|- 

Phenomenological three-nucleon forces are also employed to describe the properties of light 
nuclei with particle number up to A=8 using the Green function Monte-Carlo method [20 



However, at first sight, the situation in calculating bulk properties of finite nuclei seems to be 
different from the corresponding situation in nuclear matter. While microscopic calculations 
for nuclear matter, using the modern NN interactions, yield too much binding energy, BHF 
calculations |21| but also calculations using the coupled cluster approach p2| predict binding 
energies for finite nuclei, which are too small as compared to the experimental data. 
It is the aim of the present work to investigate this situation a bit more in detail. For that 
purpose we consider the BHF approximation paying special attention to the particle-particle 
excitations at low energies and including the effects of hole-hole scattering terms. We are 
investigating the predictions originating from various models of realistic NN interactions. 
These include local, like the Argonne V18 0, and non-local models, like the CDBonn 
interaction ||, but also the realistic meson-exchange models based on chiral perturbation 
theory, which have recently been developed by the Idaho group |23| . 

After this introduction we will describe in section 2 a technique for the calculation of bulk 
properties of nuclei which account in a consistent way for short-range and long-range cor- 
relations. Results of such calculation for the nuclei 16 and 40 Ca will be presented and 
discussed in section 3, which is followed by a short summary in section 4. 



II. CONSISTENT TREATMENT OF LONG- AND SHORT-RANGE 

CORRELATIONS 

The BHF approximation is one of the most popular approaches to account for effects of 
correlations beyond the mean-field approximation in calculating bulk properties of infinite 
nuclear matter and finite nuclei from realistic NN interactions. It is characterized by a 
solution of the Bethe-Goldstone equation, leading to the G-matrix Q, 

g,{ui) = V + V o°hq g " (1) 

and the self-consistent evaluation of the BHF single-particle energies 

*) + E (i3\G,(u = e i + e j )\ij) . (2) 

/ j<F 

In these equations V stands for the bare NN interaction, Q represents the Pauli operator, 
which restricts the propagator in the Bethe-Goldstone Eq.(|T]) to particle states with energies 
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above the Fermi energy, and i and j refer to hole-states, i.e. the eigenstates of the BHF single- 
particle hamiltonian with energies £j and Ej below the Fermi energy. The self-consistent 
definition of the starting energy uo in terms of the single-particle energies is determined by 



the Bethe-Brandow-Petschek (BBP) theorem [p4 |. Since, however, the BBP theorem can 
only be used to justify the definition of the single-particle potential for the hole states, 
the optimal definition of the single-particle energies for the particle states, which enter the 
Bethe-Goldstone Eq.([TJ) because they define the eigenvalues of H , has been discussed over 
many years in the literature. The conventional choice has been to ignore a single-particle 
potential or self-energy contributions for the particle states completely and approximate H 
by the kinetic energy only. This conventional choice is supported by the coupled cluster or 
exponential S method f25j, which using the S2 approximation essentially leads to the same 



approach, This conventional choice for the single-particle spectrum of nuclear matter is not 
very appealing as it leads to a large gap at the Fermi surface. 

Mahaux and collaborators argued that it would be more natural to choose the propagator 
according to the Green function method, i.e. define the single-particle propagator with a 
single-particle energy which includes the real part of the self-energy as a single-particle 
potential for particle and hole states. This leads to a spectrum which is continuous at the 
Fermi momentum, which provided the name "continuous choice" for this approach. This 
continuous choice leads to an enhancement of correlation effects in the medium and tends 
to predict larger binding energies for nuclear matter than the conventional choice. Inclusion 
of three-hole-line contribution [2^,p8fl indicate that the continuous choice seems to lead to 



a better convergence of the hole-line expansion and is therefore preferable. In fact, recent 
studies in nuclear matter show that the result is rather sensitive to details of the single- 



particle spectrum around the Fermi energy 29,30 



In order to investigate these effects of the low-energy particle-particle excitations, which 
should correspond to long-range correlations, also in finite nuclei, we follow the concept of a 
double partitioned Hilbert space, as it has been used before for the study of infinite nuclear 
matter (see e.g. |31| , |32| ) as well as finite nuclei (see e.g. f33"l , |34|| ). The long range correlations 
are taken into account by means of the Green function approach within a finite model space. 
This model space shall be defined in terms of shell-model configurations including oscillator 
single-particle states up certain shell. For our studies of 16 and 40 Ca we have chosen to 
include configurations up to the pf shell and sdg shell, respectively. It turned out that the 
final results are not very sensitive to this choice. This is due to the fact that the effect of 
short-range correlations, i.e. originating from configurations outside this model space are 
not ignored but taken into account by means of an effective interaction, i.e. a G-matrix 
appropriate for this model space. This effective interaction Q is determined as the solution 
of the Bethe-Goldstone equation 

G{u) = v + v — - Qm °d — g. (3) 

u ~ Vmod^oVmod 

The Pauli operator Q moc j is defined to exclude two-particle states with one of the particles 
below the Fermi level of the nucleus considered or with both nucleons in states inside the 
model space. With this G-matrix we solve the BHF Eq.(Q) using the self-consistent definition 
of the starting energy uj for hole states but also for the oscillator particle states, which are 
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inside the model space. We consider two different choices for the spectrum of the high- 
lying particle states outside the model space, defined by H in Eq.(||): The conventional 
choice, i.e. H is just the operator of the kinetic energy for the interacting particles, and 
the continuous choice, for which we add an attractive constant to the kinetic energy. This 
constant is determined in such a way, that it corresponds to the mean value of the potential 
energies of the low-lying particle states inside the model space. 

Note that this approximate BHF scheme only accounts for particle-particle correlations 
outside the model space. We will call these correlations short-range correlations in the 
discussion below. A measure for these short-range correlations is given by the depletion 
coefficient 
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or the corresponding occupation probability 

pf = l-«i. (5) 

Using the single-particle energies we define an effective interaction for the model space in 
terms of the oscillator matrix elements 

(ij \V\ kl) = ~ [(ij \Q(u = e t + ej )\ kl) + (ij \Q{u = e k + e,)| kl) ] . (6) 

The effects of long-range correlations or correlations inside the model space shall be evalu- 
ated by means of the Green function method. To determine the correlated single-particle 
Green function for a nucleon with isospin r, orbital angular momentum I and total angular 
momentum j one has to solve a Dyson equation, which in the case of a finite system with 
spherical symmetry and within a discretize model space takes the form 

g T ij(n, n'\ u) = g T ff F \n, n; u) + J2 9% (n, n"; w)AS Tli (n", n'"\ u)g Tlj (n'", n; u) , (7) 

n",ri" 

where glf HF \n, n"; uS) refers to the BHF propagator, which we assume to be diagonal in the 
radial quantum numbers n, n' 

glff F \n, n'\ u) = 6 n ,n' - . • (8) 

The correction to the BHF self-energy in terms of two particle one hole (2plh) and two hole 
one particle (2hlp) configurations inside the model space 

AE T ij(n,n',u) = S^ 1/l) (n,n',u;) + T,\j hlp \n, n' , uj) , (9) 

is defined in terms of the effective interaction V of Eq.(Q). As an example we consider the 
2plh contribution of second order in V 

v (2 P ih), , x 1 v v <nh\V\p lP2 ><p lP2 \V\n'h> 
Tlj (n ' n ' W) "2 U pi ^ >F .-(e pl + s p2 -e h ) + ir, " (10) 
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In order to evaluate the total energy of the system and expectation values of one-body 
operators one has to rewrite the single-particle Green function in the Lehmann representation 



g T ij(n,n;uj) = ^ — J - — h }^ 



(11) 



" U - LOarlj +17] ^ U) - UJ^rlj - IT] 

The single-particle density matrix is then defined in terms of the hole-spectral function 

prlj(n, n') = J2 S rlj( n i ri, UJp-rlj) (12) 

and the total energy can be evaluated from 
E 



Tlj(3n,n' 



2m 



n 



(13) 



In order to obtain the spectral function and the positions of the poles in the single-particle 
Green function, uprij, we reformulate the Dyson equation into an eigenvalue problem |2l|j3l 
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where for simplicity we have dropped the corresponding conserved quantum numbers for 
isospin, parity and angular momentum {rlj). The matrix to be diagonalized contains the 
BHF hamiltonian defined in terms of the iV BHF single-particle energies of the symmetry 
assumed within the model space, the coupling to the P different 2plh configurations and Q 
2hlp states. These couplings are expressed in terms of the matrix elements 



=< mh\V\piP2 > and A m j =< mp\V\hih 2 > 



(15) 



As long as we are still ignoring any residual interaction between the various 2pXh and 2h\p 
configurations the corresponding parts of the matrix in (|i~4t) are diagonal with elements 
in terms of single-particle energies and denoted by (Ej) for 2plh (2hlp) configurations, 
respectively. One can easily improve this approach and incorporate the effects of residual 
interactions between the 2plh configurations or 2hlp configurations. One simply has to 
modify the corresponding parts of the matrix in Eq. ( |i~4"D and replace where H.2 P ih and H.2hi P 
contain the residual interactions in the 2plh and 2hlp subspaces 




2plh 



and 




2hlp 



(16) 
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For the calculations which we are going to discuss here, we only consider the matrix elements 
for the particle-particle interaction in H.2 P ih and the hole-hole interaction for TiibXp- This 
implies that the corresponding particle-particle and hole-hole ladder diagrams are taken into 
account. A more complete treatment of the residual interaction requires the treatment of 
three-body terms and has recently been discussed e.g. by Barbieri et al. |!6] 



The eigenvalues u a of Eq. (|i4l) correspond to the poles of the single-particle Green function 
in Eq. QTTJ) and the spectral function is given in terms of the components of the eigenvectors 
by 



S h (n,n>;u; a )=X° n X° n , (17) 

for eigenvalues uo a below the Fermi energy Ep, while a corresponding equation holds for the 
spectral function S p for energies above E F . 



III. RESULTS AND DISCUSSION 

As a first example we would like to consider some results displayed in table | evaluated for 
the nucleus 16 using the CDBonn potential supplemented by the Coulomb interaction 
between protons. For the results displayed in this table the single-particle wave functions 
have been constrained to wave functions of the harmonic oscillator defined in terms of an 
harmonic oscillator constant 

(18) 



V2a ' 

with a = 0.4 fin -1 . This corresponds to an oscillator frequency of 

k, = % (19) 

of 13.27 MeV and leads to a radius for 16 0, assuming simple shell model occupancies, of 
2.65 fm, which is close to the empirical value. The first two columns of table p] refer to 
BHF calculation assuming the conventional and a continuous choice for the spectrum of the 
particle states in the Bethe-Goldstone equation ([I]). It should be recalled that we define 
the continuous spectrum in terms of the kinetic energy shifted by a constant such that the 
single-particle energies for low-lying particle states is identical to the mean value of the 
corresponding BHF single-particle energies calculated according to (0). Beside the single- 
particle energies also the occupation probabilities pf r calculated according to Eq. (0) are 
listed. Comparing these two columns one can see that the use of a continuous single-particle 
spectrum leads to more attractive single-particle energies and a binding energy, which is 
enhanced by about 2 MeV per nucleon. This enhancement of the attraction is accompanied 
by lower occupation probabilities pf r , which implies larger values for the corresponding 
depletion coefficients. This demonstrates that the lowering of the particle-state spectrum, 
going from the conventional to the continuous choice, leads to a substantial enhancement of 
correlations. 

The third and fourth column of table |, labeled BHFO and BHF(mod), refer to the model 
space approach as introduced in the preceding section. The BHFO approach identifies the 
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BHF calculation using the G-matrix denned in (|3]) with a continuous choice for the particle 
state spectrum. This means that due to the Pauli operator Q moc [ the contribution of 
low-lying particle states (for the present example those in the IsOd and lpOf shell) to the 
G-matrix are suppressed. This leads to larger occupation probabilities (compare BHFO with 
BHF (cont)) and less attractive single-particle energies. Also the binding energy per nucleon 
is reduced by about 1 MeV. 

The contribution of the low-lying particle-particle configurations are again included in the 
BHF(mod) approach by considering the corresponding (2plh) component in the definition 
of AE in Eq.(^) and solving Eq.flUp ignoring the coupling to the (2hlp) configurations but 
allowing for the residual interaction between the particle states according to Eq.([H]). In 
this way the effects of the low-lying particle-particle configurations are taken into account 
allowing for individual single-particle energies e n ij for all sub-shells with quantum numbers 
nlj and not just a replacement of these single-particle energies with the kinetic energy shifted 
by a global constant as it was done in the BHF (cont) approach discussed above. This 
more sophisticated treatment of the low-lying particle spectrum leads to some additional 
attraction in the single-particle energies ranging from 0.9 MeV in the case of the s\/2 shell 
over 1.6 MeV (P3/2) to 1.8 MeV in the case of the pi/% shell. The effect is obviously larger for 
states closer to the Fermi energy as these states are more sensitive to the details of the long- 
range correlations. The same feature can also be observed in the occupation probabilities 
Ply It should be pointed out that the results listed for the BHF (mod) approach accounts 
for depletion due to the excitation of particle state configurations inside the model space by 
means of Eq . ([12]) while the depletion due to short-range correlations leading to excitations 
outside the model space are accounted for by means of Eq.(pl). The more specific treatment 
of the low-lying particle-particle configuration reduces the spin-orbit splitting for the protons 
in the p-shell from 3.4 MeV in the case of BHF(cont) to 3.1 MeV in the BHF(mod) approach, 
both values being much smaller than the experimental one (6.3 MeV). 
For all approaches discussed so far, the single-particle strength for the hole states is just 
concentrated in one quasiparticle state. The corresponding Green functions exhibit only 
one pole below the Fermi energy. A distribution of this hole-strength is obtained if we 
also account for the 2hlp contribution in the definition of the self-energy in Eq.(|9|). This 
distribution is defined in terms of the eigenvalues u a in Eq. fll~4] ) and the strength defined in 
Eq. (|TTD . Results for the spectral distribution are displayed in Fig. [L|. In the upper panel of 
this figure, referring to the removal of protons with S1/2 quantum numbers, the strength is 
widely distributed. In the column labeled "Green" of table | we give the mean value of this 
spectral distribution defined by 

_ Y^/3n u; l3Tlj^Tlj( n T n j UJ l3Tlj) 

Tb EnPnj(n,n) 

but also the energy of the quasiparticle state (in brackets), which is defined by that eigenvalue 
^firij) which carries the largest strength. Inspecting Fig. [l| and the corresponding numbers 
in table |, one finds that the quasiparticle state for the Sx/2 state carries only rather little 
strength. 

The quasiparticle contribution is much more important for the p%/2 and pi/2 states. For these 
more weakly bound states the quasiparticle states are also those which are closest to the 
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(20) 



Fermi energy, which means that they correspond to the removal of a nucleon leading to the 
ground state or lowest excited state of the daughter nucleus of given parity (I) and angular 
momentum (J). Therefore one should consider those quasiparticle energies in comparing 
with experimental removal energies for these states (values presented in the last column of 
table D- The spin-orbit splitting resulting from these quasiparticle energies is slightly larger 
than the one derived from BHF(mod) approximation but still too small as compared with 
the experimental value. Here it should be recalled that a substantial enhancement of the 
spin-orbit splitting is obtained if the relativistic features of the Dirac Brueckner Hartree 
Fock approximation are taken into account |35fl . 

Comparing the results of the Green function approach with those obtained in the BHF(mod) 
approximation, one finds that the hole-hole terms which are included in the Green function 
approach tend to reduce the binding energy of the quasiparticle state but lead to more at- 
tractive mean values of the spectral distribution for the hole states. Since this mean value 
enters the evaluation of the total binding energy (see Eq.(|I3|)) we also get a slightly larger 
binding energy in the Green function approach as compared to the BHF(mod) approxima- 
tion. It is worth mentioning that the Green function approach with inclusion of 2hlp terms 
in the self-energy also provides non-vanishing occupation probabilities for states which are 
unoccupied in the simple shell model. Also these occupancies contribute to the total binding 
energy. 

At this stage it is useful to make a first comparison with the situation in corresponding 
calculations of nuclear matter. It has been observed also for nuclear matter that the results 
of BHF calculations are rather sensitive to the choice of the spectrum for the particle states, 
in particular for those with momenta close to the Fermi momentum p9|j30| . Using the 
precise single-particle spectrum rather than a quadratic parametrisation leads to an increase 
of the binding energy in nuclear matter around saturation density of about 1.5 MeV per 
nucleon fl30| . This is even more than the gain in binding energy from the BHF(cont) to 
the BHF (mod) approach displayed in table |. Also in nuclear matter one observes that the 
inclusion of the 2hlp terms in the self-energy lead to less attractive quasiparticle energies. 
The spread of the single-particle strength to lower energies, however, leads to additional 
binding energy. This feature, less attractive quasiparticle energies but more binding energy 
per nucleon, helps to fulfill the Hugenholtz - Van Hove theorem |37],[3£| for nuclear matter. 
The gain in binding energy due to the 2plh components in the self-energy is about 0.5 MeV 
per nucleon in nuclear matter at saturation density [|30|. Also this is a result rather similar 



to the difference between the BHF (mod) approach and the result of the complete Green 
function listed in table |. 

For a first comparison of results originating from different NN interactions, we list in table |Tl 
results of the Green function calculation for 16 using the charge- dependent Bonn interaction 
H (CD Bonn), the version A of the Idaho interaction p3j (Idaho A), which is based on 
features of chiral perturbation theory, and the Argonne V18 interaction 0. All interaction 
models yield quasiparticle energies which are more attractive than the experimental removal 
energies but a total binding energy which is slightly smaller than the empirical energy of 
7.98 MeV per nucleon. 

The prediction for the spectral strength of the p 3 / 2 and pi/ 2 quasiparticle states range be- 
tween 0.71 and 0.75. These values are larger than spectroscopic factors deduced from nucleon 
knock-out experiments, (e, e')p, which are around 0.6 [39]. The local NN interaction Argonne 
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V18 is stiffer than the non-local meson-exchange interaction models CD Bonn and Idaho 
A (see also the comparison in nuclear matter [3"0fl). This means that it predicts stronger 
correlations, as indicated by the smaller spectroscopic factors, and smaller binding energy. 
Altogether one may argue that the bulk properties of 16 O are very well reproduced. All 
interactions reproduce the total energy within 1 MeV per nucleon and predict a radius for 
the proton distribution of 2.72 fm, which is in very good agreement with the experiment. 
Here, however, one must keep in mind that the radius is to a large extent determined by the 
choice of the oscillator parameter for the model space. The choice of a = 0.4 has been made 
to obtain a radius of 2.65 fm for the uncorrelated shell-model wave function. Correlations 
within the model space lead to small enhancements, only. In order to test the sensitivity 
of the calculation on the oscillator parameter, we have performed calculations for various 
values of a. Results for the binding energy of 16 obtained in the Green function approach 
are displayed in Fig. as a function of the oscillator parameter a. These figures show that 
an energy minimum is obtained for a around 0.475 fm -1 (0.45 fm -1 , 0.44 fm -1 ) using the 
CD Bonn, the Idaho A and the Argonne V18 interaction, respectively. The corresponding 
radii of the proton distribution are 2.29 fm, 2.41 fm and 2.48 fm, all of them to small as 
compared to the experimental value. Therefore the situation for calculating bulk properties 
of 16 is similar to the attempts of evaluating the saturation point of nuclear matter. The 
calculation tend to predict a radius which is to small or a density which is too large. It is 
worth noting, however, that the disagreement for nuclear matter is larger. 
As a second example we also consider the nucleus 40 Ca. In this case we consider a model 
space, which is defined in terms of harmonic oscillator wave functions including all shells up 
to the 2sld0g shell. For a first comparison we fix the oscillator constant to a = 0.35 fm -1 . 
This corresponds to simple shell model prediction for the radius of 3.5 fm, in reasonable 
agreement with the empirical value of 3.48 fm. 

Examples for the spectral distribution are displayed in Fig. |3] (assuming the CD Bonn 
interaction) and numerical results obtained in the Green function approach are listed in table 
PI . Since we have two quasiparticle states for the sin channel in the BHF approximation 
(see upper panel of Fig. ||), we also list two energies and occupation probabilities for this 
partial wave in table (TTT] . Also for this nucleus one obtains a broad distribution of the spectral 
strength. A well defined quasiparticle peak only shows up for the states in the IsOd shell. 
Also in this case, the predictions for the removal energy (absolute value of the quasiparticle 
energy for (i 3 / 2 , are larger than the experimental value of 8.3 MeV. Note, however, that the 
inclusion of 2hlp contributions in the Green function approach reduces the discrepancy by 
more than 3 MeV as compared to the BHF (mod) approximation. On the other hand, the 
calculated binding energies are in good agreement with the experimental value of 8.55 MeV 
per nucleon. So again, the calculation of bulk properties of 40 Ca is rather successful, if one 
fixes the radius with the appropriate choice of the oscillator constant a, which defines the 
basis of the model space. 

If one releases this constraint and considers various values for a one obtains results as 
displayed in Figs. [| and ^. Also in this case we observe that the minima occur for oscillator 
parameters which are larger than a = 0.35 fm -1 , which means that the corresponding radii 
are smaller than the empirical one. This is true for the various approximations (see Fig. |J) 
but also for the various interactions (see Fig. 0). 

The comparison of the various approximation schemes in Fig. H confirms the conclusions, 
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which have been given above for the case of 16 0. In the case of Ca, however, all interactions 
predict too much binding energy and a radius, which in the case of CD Bonn interaction 
is about 25 percent smaller than the experimental value. This corresponds to an average 
density, which is too large by about a factor of two, a situation which is approaching the 
situation of nuclear matter. Also in this case, the stiffer interaction Argonne V18 yields less 
energy and smaller densities than the softer interactions (Idaho A and CD Bonn). 



IV. CONCLUSIONS 

Starting from the Brueckner-Hartree-Fock (BHF) approach various approximation schemes 
have been investigated to derive bulk properties of finite nuclei from realistic NN interactions. 
It is observed that the results of BHF calculations, for finite nuclei as well as infinite nuclear 



matter |[30|| , are rather sensitive to the spectrum of particle states, in particular those with 
energies close to the Fermi energy. Therefore a technique has been applied, to separate a 
detailed description of the low-energy excitations corresponding to long range correlations 
from the treatment of short-range correlations. While the effects of short-range correlations 
are taken into account by means of the G-matrix approximation of the BHF scheme, the 
long-range correlations are considered within the framework of the self-consistent evaluation 
of single-particle Green function. This approach includes the effects of particle-particle 
correlations but also corresponding hole-hole scattering terms. 

If the basis of the model space is constrained to obtain the empirical value for the radius, one 
is able to reproduce the binding energy of nuclei (we consider 16 and 40 Ca) within 0.5 MeV 
per nucleon, which is very satisfactory for a calculation based on realistic NN interaction. 
This success is due to the attraction which is obtained from the careful treatment of the 
low-lying particle-particle excitations but also due to the inclusion of two-hole one-particle 
configurations in the definition of the self-energy for the single-particle propagator. These 
terms yield a distribution of the single-particle strength and a shift of the quasiparticle 
energy. The repulsive shift of the quasiparticle energy for states close to the Fermi level 
improves the agreement with the experimental removal energies. The distribution of the 
strength, on the other hand, allows a gain in binding energy although the quasiparticle 
energies are less attractive. This helps to improve the fulfillment of the Hugenholtz - Van 
Hove theorem. The calculated spectroscopic factors for the quasiparticle states (around 0.7) 
are slightly above the empirical values (around 0.6) derived from (e, e')p experiments. 
If the constraining condition for the radius is released, the minima for the calculated energies 
lead to radii which are too small and densities which are too large. The calculated binding 
energies tend to be larger than the empirical values. This is more pronounced for the nucleus 
40 Ca than for 16 0. This deficiency could be cured by considering the repulsive three nucleon 
forces f|-12], which are adjusted to describe the saturation point of infinite nuclear matter 
also for the calculation of finite nuclei. 
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BHF (conv) 


BHF (cont) 


BHFO 


BHF (mod) 


Green 


Exp 


Si./ A 


-38.19 


-42.78 


-40.74 


-43.72 


-44.00 (-47.01) 


-40 ± 8 


£ P3/2 


-18.14 


-22.40 


-20.40 


-23.99 


-24.29 (-20.68) 


-18.45 


£ Pl/2 


-14.50 


-19.02 


-16.86 


-20.86 


-21.26 (-17.44) 


-12.13 


Psl/2 


0.928 


0.854 


0.904 


0.871 


0.867 (0.276) 




Pp3/2 


0.926 


0.855 


0.898 


0.829 


0.823 (0.730) 




Ppl/2 


0.916 


0.857 


0.900 


0.818 


0.802 (0.725) 




E/A 


-4.61 


-6.67 


-5.66 


-7.57 


-7.78 


-7.98 



TABLE I. Single-particle energies (ej) and occupation probabilities (pi) for protons in 16 and 
the total energy per nucleon (E/A) as calculated from various approximations (see discussion in 
the text) are compared to experimental data. The calculations have been performed using the 
CD Bonn interaction and considering a model space defined in terms of oscillator functions with 
parameter a = 0.4. The numbers in brackets refer to energy and occupation probability of the 
dominant quasiparticle contribution in the Green function approach. All energies are given in 
MeV. 
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CD Bonn 


Idaho A 


Arg. V18 


e sl/2 


-44.00 (-47.01) 


-44.00 (-47.43) 


-42.38 (-45.63) 


£ P'i/2 


-24.29 (-20.68) 


-24.00 (-20.64) 


-23.12 (-19.77) 


£ pi/2 


-21.26 (-17.44) 


-20.87 (-17.28) 


-20.26 (-16.62) 


Pal/2 


0.867 (0.276) 


0.873 (0.284) 


0.830 (0.265) 


Pp3/2 


0.823 (0.730) 


0.833 (0.747) 


0.806 (0.723) 


Ppl/2 


0.802 (0.725) 


0.821 (0.741) 


0.794 (0.715) 


E/A 


-7.78 


-7.65 


-7.15 



TABLE II. Results for single-particle energies, occupation probabilities and total energy per 



nucleon of O calculated from the Green function approach using the CD Bonn, the Idaho A and 
the Argonne V18 interaction model. Further details see table |. 





CD Bonn 


Idaho A 


Arg. V18 


e sl/2 


-38.30 (-57.26) 


-38.81 (-58.44) 


-37.92 (-55.53) 




(-17.63) 


(-17.94) 


(-16.84) 


£ p 3 /2 


-37.59 (-23.22) 


-38.06 (-23.76) 


-36.37 (-22.52) 


£ Pl/2 


-35.17 (-22.16) 


-35.60 (-22.59) 


-34.12 (-21.42) 


£ d 5 /2 


-22.27 (-18.35) 


-22.42 (-18.70) 


-21.33 (-17.61) 


£ d 3 /2 


-18.88 (-15.08) 


-18.96 (-15.34) 


-18.13 (-14.46) 


Psl/2 


1.658 (0.160) 


1.674 (0.177) 


1.592 (0.142) 




(0.697) 


(0.715) 


(0.682) 


Pp3/2 


0.834 (0.123) 


0.844 (0.141) 


0.800 (0.151) 


Ppl/2 


0.826 (0.301) 


0.838 (0.309) 


0.794 (0.292) 


Pdb/2 


0.813 (0.691) 


0.824 (0.709) 


0.795 (0.679) 


Pd3/2 


0.796 (0.701) 


0.808 (0.717) 


0.780 (0.688) 


E/A 


-8.77 


-8.91 


-8.22 



TABLE III. Results for single-particle energies, occupation probabilities and total energy per 



nucleon of Ca calculated from the Green function approach using the CD Bonn, the Idaho A and 
the Argonne V18 interaction model. Further details see table |. 
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Energy [MeV] 

FIG. 1. The spectral function for proton hole strength in the S1/2 (upper panel) andj>3/2 channel. 
The results are calculated for 16 using the CD Bonn interaction. The distributions of the complete 
Green function approach are obtained by folding the discrete distribution with Gaussian functions 
assuming a width of 1 MeV. Also given are the positions of the single-particle states evaluated in 
the BHF(mod) approximation (straight lines). 
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_9 1 1 1 1 1 

0.4 0.45 0.5 

Oscillator parameter a [fm ] 

FIG. 2. Binding energy per nucleon calculated in the Green function approach as a function of 
the oscillator parameter a (see Eq.( |l8|), which is used to define the basis of the model space). 
Results are displayed for various NN interactions 
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-60 -40 
Energy [MeV] 



FIG. 3. The spectral function for proton hole strength in the S]/ 2 (upper panel) and pi/ 2 (middle) 
and d 3 / 2 channel. The results are calculated for °Ca using the CD Bonn interaction. Further 
details see Fig. [fl. 
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■10 - 



40 Ca 

Argonne 



-- BHFO 

■ ■ ■ BHF(mod) 

— Green 



0.36 



0.4 



0.44 



0.48 



Oscillator parameter a [fm ! ] 



FIG. 4. Binding energy per nucleon for 40 Ca calculated in the Green function approach as a 
function of the oscillator parameter a. Results are displayed for various approximation schemes 
(for nomenclature see table |) 
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-12 1 — 1 1 1 1 1 1 1 

0.36 0.4 0.44 0.48 

Oscillator parameter a [fm *] 

FIG. 5. Binding energy per nucleon calculated in the Green function approach as a function of 
the oscillator parameter a. Results are displayed for various NN interactions 
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